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ABSTRACT 

We  consider  a  viscoelastic  fluid  filling  a  bounded  domain  in  R^  under 
the  influence  of  a  small  body  force.  The  fluid  is  described  by  certain 
differential  constitutive  equations.  We  use  an  iterative  method  to  prove  the 
existence  of  steady  flows. 
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EXISTENCE  OF  SLOW  STEADY  FLOWS  OF  VISCOELASTIC  FLUIDS 
WITH  DIFFERENTIAL  CONSTITUTIVE  EQUATIONS 

Michael  Renardy 

1.  INTRODUCTION 

Although  existence  theorems  for  low  Reynolds  number  steady  flows 
of  Newtonian  fluids  are  well  known  (see  e.g.  [5]),  no  such  theorems  have 
been  established  for  viscoelastic  fluids.  There  are  formal  perturbation 
expansions  which  are  used  for  slow  flows  (see  e.g.  [10]),  but  the 
justification  of  these  expansions  leads  to  difficult  problems  in  singular 
perturbation  theory,  which  have  not  been  solved.  Niggemann  [8]  has 
given  a  convergence  proof  for  expansions  of  this  nature  in  a  one-dimensional 
model  problem,  which  has  certain  features  in  common  with  equations  in 
viscoelasticity. 

In  this  paper,  we  prove  the  existence  of  slow  steady  flows  of  certain 
viscoelastic  fluids  by  using  an  iterative  method.  The  basic  idea  is  very 
similar  to  existence  proofs  for  initial  value  problems  in  hyperbolic  partial 
differential  equations.  We  first  show  that  all  iterates  are  bounded  and 
small  in  a  certain  norm,  and  we  then  show  that  the  iteration  converges 
in  a  weaker  norm.  The  iteration  we  use  is  similar  to  procedures  employed 
in  numerical  calculations  (for  a  review,  see  [1]),  and  the  ideas  used  here 
should  therefore  be  useful  in  proving  convergence  of  numerical  schemes. 

We  shall  also  use  our  results  to  justify  the  formal  perturbation  methods 
as  asymptotic  expansions,  but  we  do  not  prove  their  convergence. 


We  study  steady  flows  of  a  viscoelastic  liquid  in  a  bounded  donain 
ficj3.  The  fluid  satisfied  the  no-slip  condition  on  the  wall  and  moves 
under  the  influence  of  a  given  body  force.  (Problems  arising  in  applications 
usually  have  inflow  boundaries,  which  require  additional  boundary  conditions 
Such  problems  are  more  difficult  than  the  one  studied  here) .  The  equations 
for  steady  flow  are  as  follows 


(1.1) 


p(u*V)u  "  div J  ♦  Vp-f  «  0 

divii  ■  0 

u  »  0 


in  (ICR3 
on  . 


Here  u  *  (Uj.u^u^)  is  the  velocity  vector,  p  the  pressure  and  T 
the  extra  stress,  f  is  the  given  body  force  and  p  is  the  density. 
Throughout  the  paper,  12  is  assumed  to  be  a  bounded  domain  with  a 
smooth  (for  simplicity,  say  C“-)  boundary. 

The  extra  stress  is  related  to  the  velocity  field  by  a  constitutive 
equation.  Here  we  deal  with  differential  constitutive  equations.  For 
sinqplicity,  we  adopt  a  particular  constitutive  law,  which  exemplifies 
the  typical  structure.  This  model,  the  "rubberlike  liquid"  [3],  [7], 
is  given  by 


N 

I  •  2"oS  ■»  I  Xk  . 

k«l 

(“•'Oik  -  TOIk  -  Ik(*i)T  »  *kXk  ■  Vk8  • 


(1.2) 


and 


where  Vu  is  the  velocity  gradient  :  (Vu). 
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D  ■  |-(Vu  ♦  (Vu)T).  N  is  an  arbitrary  positive  integer  and  are 

positive  constants  (nQ  may  be  zero). 


The  essential  difficulty  for  the  analysis  arises  from  the  teras 

T 

(u'V)!^.  The  particular  form  of  the  terms  (Vu)T^  and  T^(7u)  is 
unimportant ,  and  we  could  replace  them  by  other  nonlinear  combinations 
of  (Vu)  and  T^.  Our  analysis  can  thus  be  extended  to  fluids  with 
other  differential  constitutive  equations,  such  as  those  of  Oldroyd  (9] 
Leonov  [6]  and  Giesekus  [2] . 


2.  THE  CASE  OF  THE  UPPER  CONVECTED  MAXWELL  aUID 


In  this  case,  we  set  i)g  •  0  and  N  ■  l  ,  i.e.  the  constitutive 


equation  takes  the  fora 


(2.1) 


(u*7)|  -  (7u)T  -  T(Vu) 1  ♦  XJ  -  2nXg 


By  applying  the  divergence  operator,  we  obtain 


(2.2) 


(u*7)divT  -  (Vu)div  T  ♦  AdivT  «  J  :  3*u  ♦  qAAu  . 


Here  we  use  the  notation  T  :  3*  «  Y  T..  * — 5 —  .  Next,  we  substitute 

-  j>k  3xjaxk 

div  T  »  p(u*V)u  ♦  Vp  -  f  from  (1.1),  and  obtain 


(2.3) 


7[(u*7)p  +  Xp]  -  [(7u)  ♦  (Vu)  ]7p-  (u*7)f 
♦  (7u)f—  Xf  «  T  :  32u+nXAu  -  p(u*7)(u*7)u 


*  p(7u) (u*7)u- Ap(u*7)u 


In  the  following,  we  regard  (2.3),  with  the  condition  div  u-0 
and  the  no>slip  condition,  as  a  perturbation  of  the  Stokes  problem.  This 
equation  contains  a  "modified  pressure"  q  *  (u*7)p  ♦  Xp  .  Solutions  are 
found  by  the  following  iteration  scheme. 


(2.4) 


(2.S) 


0  _0 


uo-g,p-q«0,  I-g 


t”  :  3zun+1  ♦  nXAun+1  -  p(un,7) (un,7)un+1  -  7qf 


n.m  /. 


-  - [7un  ♦  (7un)T]7pn—  (un*V)f  ♦  (7u”)f — Xf 


-  P(7un) (un»7)un  ♦  Ap(un *7)un  , 


div  un+1  ■  0  ,  un+1 *  0  on  3fl  ,  ///  q”*1  -  0  , 

n 


I 


vV.TV 

*  *  «  * 
••V'Vv 


;•  Vv 


t.V/.'.N 


w  ■ 

\V 

Wsi1 

t 

Ws 
.  •  ^  ^ 


h*  '•*  v", 


v.vv 


(2.7) 


(un+1*V)Tn+1  -  (7un+1)Tn+1-  Tn+1(Vun*Y 


♦  ATn+1  *  nA(7un+1  ♦  (7un+1)T)  . 


We  denote  by  Hs(fl)  the  usual  Sobolev  spaces  and  by  IMIS  the 
norm  in  Hs(fl)  .  The  following  lemma  is  immediate  from  the  invertibility 
of  the  Stokes  operator  (SJ  and  elementary  perturbation  theory  [4J. 

LEMMA  2.1  :  Let  s  be  an  integer  >  1  .  Then  there  are  positive  constants 

and  e2  such  that  the  following  holds:  If  HX*l|ls+1  -  ^  ora? 

||un|ls+1  <  e2  ,  then  equation  (2.S)  has  a  unique  solution  un+1  ,  q11*1  . 

This  solution  obeys  an  estimate  of  the  form 


hu"*1iis,2.  ii<r1ns>1  s  c1dip”iis.1iis"iis.2. 


Bflls.||u"llJ.2||>!nlls*llyn|l'.i] 


The  following  lemma  concerns  the  solvability  of  equation  (2.6). 

LEMMA  2.2:  There  is  some  >  0  such  that,  for  ||un+*||s+2  -  e3  •* 

div  un+1  *  0  ,  un+1|  *  0  ,  the  unique  solution  pn+1  of  (2.6)  satisfies 


an  estimate  of  the  form 


l!pn+1lts+1  s  c2||,"*1||J.1 


Sketch  of  the  proof:  Equation  (2.6)  is  easily  solved  by  the  method  of 
characteristics,  and  existence  and  uniqueness  follow  immediately.  Moreover, 
the  operator  (u  *V)  is  skew-adjoint  in  L ^(fl)  ,  whence 

„  n+1,,  1  n*l„ 

IIP  H0  5  7  IW  M0  * 

A 


vow 


*  V ’•/  An-’  •’ 


Estimates  for  derivatives  are  obtained  by  differentiating  the 
equation.  (Such  estimates  are  formal  and  the  calculations  involve  derivatives 
not  a  priori  known  to  exist.  However,  it  is  easy  to  see  that  p  is 
smooth  if  u  and  q  are,  and  we  can  thus  construct  approximating  sequences 
satisfying  uniform  estimates.) 

Since  equation  (2.7)  can  be  regarded  as  a  perturbation  of  (2.6), 

we  have 

LEMMA  2.3:  There  is  a  constant  c ^  such  that ,  for  ||un**||s+2  ^ 
unique  solution  of  (2.7)  satisfies  an  estimate  of  the  form 

»r\.i  *  • 

By  combining  the  estimates  contained  in  lemmas  2.1  -  2.3,  it  can 
be  shown  that  all  iterates  remain  bounded  if  f  is  small. 

LEMMA  2.4:  If  ||f|ls+1  is  sufficiently  small ,  then  ||un||s+2  ,  l|pnHs+1  »• 
liqn„s+l  and  ll^llj+i  haue  (small)  bounds  independent  of  n  . 

Next  we  show  that  the  iteration  generates  a  convergent  sequence  in 
a  weaker  norm. 

LEMMA  2.5:  Let  llflls+1  be  sufficiently  small.  Then  there  is  a  constant 
Y  <  1  such  that 

Hu"*'1  -  u”|jn  ♦  Iq"*1  -  qn»5  *  Y(|un  -  un'l||  Jtl  ♦  llq“-  I*'1!!,!  • 

Sketch  of  the  proof:  From  (2.6),  we  obtain 

,  n  n  ,  n-1  _  n-1  , ,  n  n-l.  n  n-1 

(u  *7)p  -  (y  *7)p  ♦  Mp  -P  )  ■  q  -  q 


This  is  equivalent  to 


((un-un~Vv)pn  ♦  (un'1*7)(pn-pn"1)  ♦  X(pa-p"’1) 


n  n-1 

q  -  q 


From  this  and  the  bounds  already  established  by  1 
conclude  that,  for  some  constant  ,  we  have 

llpn-pn‘1|ls  SC4[||qn-qn-1|ls  ♦  Il«n-Un-I|ls+Il  • 


2.4,  we 


Similarly,  we  find  from  (2.7) 

l|Tn-Tn"1| 


s  2CSH“n-“n'1,,s4l  ‘ 


Next,  we  subtract  (2.5)n  and  (2.5)n_j  .  By  using  lemma  2.4  and 
the  already  established  estimates  for  pn  -  pn  *  and  J11-  ,Jn  *  ,  we  easily 
obtain  the  lemma.  We  omit  the  details  of  the  calculation. 

Thus  we  have  proved 

THEOREM  2.6:  Let  s  be  an  integer  >  1  and  let  ll?llJ+1  be  sufficiently 
small.  Then  there  exists  a  solution  u  €  HS+*  ,  p  €  H*+*  ,  T  €  HS+*  for 
equations  (1.1),  (2.1),  obtainable  by  the  iteration  procedure  (2.4)-  (2.7). 

Remark: 

Let  us  replace  f  by  ef  .  If  s  is  chosen  large  enough,  then  we 
can,  by  following  similar  procedures  as  above,  obtain  estimates  for 
difference  quotients  of  the  solution  with  respect  to  c  .  This  shows  that 
the  solution  depends  smoothly  on  e  and  therefore  establishes  the 
asymptotic  validity  of  Rivlin-Ericksen  expansions.  In  problems  with 
inflow  boundaries,  we  should  expect  the  situation  to  be  quite  different. 
Rivlin-Ericksen  expansions  are  uniquely  determined  by  prescribing  velocities 
alone  on  the  boundary.  However,  these  boundary  conditions  are  clearly  not 
enough  to  uniquely  define  flows  of  a  Maxwell  fluid. 
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We  apply  the  divergence  operator  to  (1.2),  and  obtain  as  before 

(3.1)  (u*7)div  Tk-  7u  div  |k*  div  |k 

2 

*  nkXkAu  ♦  Tk  :  d  u  . 

We  write  this  as 

(3.2)  dxv  Tk  «  ((u*V)  ♦  Xk)_1[nkXkAu*Tk  :  32u+7u  div  |k)  . 

By  inserting  this  into  (1.1),  we  find 

(3.3)  p(u*7)u  «  nnAu  ♦  I  ((u*7)  ♦  X.  )"* 

'  "  u  "  k  '  * 

[nkXjAu+Tk  :  32u  ♦  7u  div  Tfc]  -  7p+  f  . 

The  cases  nQ  t  0  and  nQ  M  0  are  treated  in  different  ways.  We  begin 
with  n0  J4  0  •  I"  this  case,  we  set 

L[u]  -  n0  ♦  X  nkXk((u«V)  ♦  Xk)_1  . 

LEMMA  3.1:  The  operator  L[u]  is  a  bisection  L2(fl)  L2(Q)  . 

Proof:  Since  (u*7)  is  skew-adjoint,  we  have  (((u*7)  ♦  Xk)“*x,x)  >  0 

x  €  l2(J1)  .  Hence  (L(u]x,x)  >  Hq(*,x)  ,  and  the  invertibility  follows 
from  the  Lax-Mil gram  theorem. 

By  differentiating  the  equation  L[u]x*y  ,  it  is  not  difficult 
to  show  that,  if  u  and  its  derivatives  are  small,  then  L[u]  also  maps 
higher  Sobolev  spaces  bijectively  into  themselves. 

I 

We  set  q  *  l[u)"lp  and  apply  the  operator  l[u]~l  to  (3.3). 


We  thus  obtain 


(3.4) 


pL[u]  (u-V)u  =  A u-  Vq ♦  L[u]_1  {  £  ((u-V)  +  X.  )'A 
‘  ’  ‘  '  k 

(Tk  :  32u  +  7u  div  TjJ)  +  L[u)_1f 
♦  L[u) _1  l  nkXk((u-V)  +  Xk)"1(Vu)TV[((u-V)^Xk)"1q]  . 


We  use  the  following  iteration  scheme 


(3.S) 

u°  *  0  ,  q°  =  0  ,  Tk  *  0 

(3.6) 

.  n+1  „  n+1  .  r  n.-l,  n  n 

Au  -  7q  =  pL[u  ]  (u  *v)u 

i r  n.-l 
-  L[u  ] 

l  ((un*V)  +  Xk)_1[TjJ:  32un+ VundivTjJ] 
k 

. r  n.-l 
-  L[u  ] 

f  -  L[un]-1 1  nkXk((u"-V)  +  Xk)_1(7un)T 

v[((un-v)  ♦  \y\n) 

n+1 
div  u 

-  0  ,  un'“1|  =  0  ,  ///  qn+1  =  0  . 

3ft  ft 

(3.7) 

,  n+1 
(u 

•V)ljJ+1  -  (Vu)n+1Tj+1  -  Tj+1(Vun+1)T 

+ 

XkS+1  =  \Xk[Vun+1+  (7un+1)T)  . 

We  can  now  proceed  in  precisely  the  same  manner  as  in  section  2  to 
establish  convergence  of  the  iteration  scheme. 

If  Hq  =  0  ,  the  operator  L[u]  as  defined  above  is  not  coercive 
and  lemma  3.1  does  not  hold.  In  this  case,  we  adopt  a  different  procedure 
Let  X  be  any  positive  real  number.  We  apply  the  operator  (u»7)  +  X  to 


(3.3)  and  obtain 


(3.3)  p((u*7)  ♦  X)  (u*7)u  »  ((u*7)  ♦  X){  £  ((u*7)  ♦  X.  )‘1 

"  "  '  k 


[f^Xk^u  ♦  Tk  :  3Zu  +  Vu  div  Tk]  -  7p  ♦  £  } 


We  now  set  L[u]  ■  ((u*7)  ♦  X)  £  nk*k  ((“‘^)  ♦  ‘ 

*  I  nkxkU*  (X-Xk)((u-7)>Xk)'1}  .  . 

2  -1 

L(u]  is  a  coercive  operator  in  L  (ft)  and  L[u]  exists.  We 

-1  _  i  ~ 

can  now  apply  the  operator  L[u]  to  (3.8),  define  q=  L[u]  ((u»7)  ♦  X)p 
and  set  up  an  iteration  scheme  in  an  analogous  fashion  as  before. 
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20.  ABSTRACT  ( Continue  on  rover  me  elde  II  necessary  md  Identity  by  block  number) 

We  consider  a  viscoelastic  fluid  filling  a  bounded  domain  in  ]RJ  under 
the  influence  of  a  small  body  force.  The  fluid  is  described  by  certain 
differential  constitutive  equations.  We  use  an  iterative  method  to  prove  the 
existence  of  steady  flows. 
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